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ABSTRACT 

We reconstruct the history of the cosmic star formation in the universe by means of detailed chemical 
evolution models for galaxies of different morphological types. We consider a picture of coeval, non- 
interacting evolving galaxies where ellipticals experience intense and rapid starbursts within the first Gyr 
after their formation, and spirals and irregulars continue to form stars at lower rates up to the present 
time. Such models allow one to follow in detail the evolution of the metallicity of the gas out of which 
the stars are formed. We normalize the galaxy population to the B band luminosity function observed 
in the local Universe and study the redshift evolution of the luminosity densities in the B, U, I and 
K bands calculating galaxy colors and evolutionary corrections by means of a detailed synthetic stellar 
population model. Our predictions indicate that the decline of the galaxy luminosity density between 
redshift 1 and observed in the U, B and I bands is caused mainly by star-forming spiral galaxies which 
slowly exhaust their gas reservoirs. Elliptical galaxies have dominated the total luminosity density in all 
optical bands at early epochs, when all their stars formed by means of rapid and very intense star-bursts. 
Irregular galaxies bring a negligible contribution to the total luminosity density in any band at any time. 
We study the cosmic missing metal crisis and conclude that it could be even more serious than what has 
been assessed by previous authors, if the bulk of the metals were produced in dust-obscured starbursts 
associated with the early spheroid formation. The most plausible site for the missing metals could be 
the warm gas in galaxy groups and proto-clusters, in which the metals could have been ejected through 
galactic winds following intense starbursts. Finally, we predict the evolution of the cosmic star formation 
and supernova II and la rates and obtain the best fit to the observations assuming a Salpeter IMF. All 
our results indicate that the bulk of the stellar mass in most galaxies was already in place at early epochs 
and we predict the existence of a peak in the global SFR at high redshift due to ellipticals. 

Subject headings: galaxies: evolution, galaxies: fundamental parameters, galaxies: photometry, 
galaxies: high-rcdshift 



1. INTRODUCTION 

The study of the evolution of the galaxy luminosity den- 
sity is fundamental to understand how galactic structures 
formed and evolved in the universe. The light emitted by 
stars of various masses at various wavelengths can provide 
different indications concerning how the cosmic star for- 
mation history proceeded in the past and concerning the 
fraction of baryons locked up in stars and gas in the local 
universe. The short wavelength light, i.e. that emitted 
in the rest-frame U and B bands, is mainly radiated by 
short-lived massive stars, thus it can be a direct tracer of 
star formation. On the other hand, the near-infrared light 
is primarily emitted by low mass, long living stars which 
contribute to the bulk of the total stellar mass; hence it 
traces the mass ditribution of galaxies. Several attempts 
to model the cosmic history of star formation and the evo- 
lution of the luminosity density have been performed, gen- 
erally following different routes. One way to reconstruct 
the evolution of galaxies is often referred to as the "tra- 
ditional" scheme, where galaxy densities are normalized 
according to the optical and IR luminosity functions ob- 
served at 2 = 0. The adoption of a star formation history 
and a cosmological model allows one to follow the redshift 
evolution of the luminosity density as well as to recon- 
struct other properties, such as number counts and the 
color distributions (Totani, Yoshii, Sato 1997; Pozzetti et 
al. 1998; Tan, Silk & Balland 1999; Jimenez & Kash- 



Hnsky 1999; Totani & Takeuchi 2002). A different ap- 
proach is followed by Madau, Pozzetti & Dickinson (1998), 
who start from the observed time-dependent star forma- 
tion rate (SFR) per unit comoving volume and an initial 
mass function (IMF), and then model the luminosity den- 
sity by means of a photometric code. Fall, Chariot & Pei 
(1996) compute the evolution of the total cosmic emissiv- 
ity and of the background intensity on the basis of quasar 
absorption-line studies, i.e. HI and metallicity data ob- 
served in Damped Lyman- alpha (DLA) systems, combined 
with cosmic chemical evolution (Pei & Fall 1995) and stel- 
lar population models (see also Sadat, Guiderdoni & Silk 
2001 and Bosclli ct al. 2001 for similar approaches). Other 
groups compute the cosmic star formation history using 
large-scale hydrodynamical simulations (Nagamine, Cen, 
Ostriker 2001, Ascasibar et al. 2002) or semi-analytical 
models of galaxy formation (Baugh et al. 1998, Cole et 
al. 2000, Somerville, Primack, Faber 2001), both methods 
based on the cold dark matter paradigm. An interesting 
alternative is the approach by Rowan-Robinson (2001), 
where the star formation rate is parametrized in order to 
investigate the effects that the choice of these parameters 
has on the far-infrared and submm counts and background 
radiation. 

However, most of the approaches described above cannot 
predict the roles of the different galactic morphological 
types at various cosmic epochs, and which type of galaxy 
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determines the strong evolution observed in the luminosity 
density between z = and z = 1. 

In this paper, we calculate the evolution of the galaxy lu- 
minosity density and cosmic star formation rate by means 
of detailed chemical evolution models for galaxies of dif- 
ferent morphological types, i. e. ellipticals, spirals and 
irregulars, which successfully reproduce the local proper- 
ties of such galaxies. We match our chemical evolution 
models with a, photometric code which allows us to com- 
pute the galaxy spectra and magnitudes, when the IMF is 
fixed. We normalize the galaxy fractions according to the 
B-band luminosity function (LF) as measured by Marzke 
et al. (1998), and we compute the LF in other bands 
as a function of redshift on the basis of the photometric 
evolutive corrections and colors predicted by our models. 
Our detailed study allows us to provide an answer to sev- 
eral questions concerning the cosmological evolution of the 
baryons in the universe in the form of stars and metals. In 
particular, we study the cosmic metal production and the 
related missing metal problem. This can indicate how the 
chemical enrichment of the inter-galactic medium (IGM) 
has occurred and what might be the nature of the objects 
observed in the high-redshift universe, such as Lyman- 
break (LB) and hyperluminous infrared galaxies. Finally, 
we focus on the cosmic evolution of type la and II super- 
nova (SN) rates, and gain further independent probes of 
how galactic structures have evolved in the universe since 
their formation. 

The paper is organized as follows: in section 2 we review 
the current observational status concerning the luminosity 
density and the star formation measurements in low- and 
high-redshift galaxies, in section 3 we describe the theory 
at the basis of our chemo-spectrophotometric models, in 
section 4 we present our results and in section 5 we draw 
the conclusions. 

2. OBSERVATIONS OF HIGH-REDSHIFT GALAXIES: 
INTEGRATED LUMINOSITY DENSITY AND COSMIC 
STAR FORMATION RATE 

2.1. Determining the galaxy luminosity function 

The total luminosity density (hereinafter LD) in a given 
band is the integrated light radiated per unit volume from 
the entire galaxy population. The distribution of absolute 
magnitudes for galaxies of any specified Hubble type is 
represented by the himinosity function (LF, Efstathiou et 
al. 1988, Binggeli et al. 1988). The luminosity function 
is often parametrized according to the form defined by 
Schechter (1976): 

$(L) dL/L* = (L/i*)-" exp{-L/L*) dL/L* (1) 

where is a normalization constant related to the num- 
ber of luminous galaxies per unit volume, L* is a charac- 
teristic luminosity and a is associated to the slope of the 
function, which accounts for the percentage of faint sys- 
tems. The LD stems from the integral over all magnitudes 
of the observed luminosity function: 

PL = J ^L){L/L*)dL (2) 

Different determinations of the local field LF (Loveday et 



al. 1992, Marzke et al. 1994b, Elhs ct al. 1996, Marzke 
et al. 1998, Cross et al. 2001) show several discrepancies 
(Wright 2001), in the sense that all the three parameters 
are not fully constrained by observations (Ellis 1997). Re- 
cent progress has being made in the determination of the 
LF as a function of morphology (Marzke et al. 1994a, 
Marzke et al. 1998, Brinchmann et al. 1998, Kochanek et 
al. 2001) and of spectral type (Heyl et al. 1997, Madgwick 
et al. 2002) or color (Blanton et al. 2001), which is impor- 
tant in determining the relative contributions of different 
galaxian types to the total luminosity density in the lo- 
cal universe. The contributions of different morphological 
types is a function of the band in which the LF is esti- 
mated, with early-type galaxies contributing a substantial 
fraction of the total emissivity at long wavelenths, namely 
in the I and K bands where old massive galaxies domi- 
nate (Kochanek et al. 2001), whereas late star-forming 
spirals arc the major contributors in the B and U bands 
(Marzke et al. 1998, Fukugita et al. 1998). Notwithstand- 
ing the uncertainties in the LF parameters, there is an 
overall concordance among various observations regarding 
the redshift evolution of the galaxy LD. This overall agree- 
ment concerns, in particular, the < z < 1 redshift range, 
where the LD is observed to rise sharply, although the pre- 
cise steepness is still under debate (Cowie et al. 1999, Lilly 
et al. 2002). On the other hand, the high redshift trend 
is rather uncertain, since surface-brightness dimming ef- 
fects could be significant (Lanzetta ct al. 2002) and since 
the bulk of the available data has been performed in the 
rest-frame UV, which can be seriously affected by dust ob- 
scuration effects. Dust tends to absorb the UV light emit- 
ted by young stars and to re-radiate it in the IR-submm 
bands. Moreover, the extent to which the presence of dust 
contaminates the data is rather uncertain, since galaxies 
with very different bolometric luminosities can have very 
similar UV luminosities (Adclbcrgcr 2001). Some authors 
(Connolly et al. 1997, Madau et al. 1998) claim a broad 
peak in the UV LD located at z > 1, followed by a fall 
at higher redshifts. Other authors observe a rather con- 
stant behaviour at z > 1 (Sawicki et al. 1997, Pascarelle 
ct al. 1998), which extends to an imccrtain epoch when 
the first galactic structures formed. In particular, accord- 
ing to the recent Hubble Deep Field data, the constancy 
of the galaxy LD is observed out to z 6, as reviewed in 
Thompson (2002). 

2.2. The observed star formation rate density 

The determination of the SFR density is related to the 
measures of star formation in galaxies, which can be per- 
formed by means of various emission processes (Madau 
1997, Kennicutt 1998a, Schacrcr 1999). AU the different 
measures require the assumption of a universal mass func- 
tion, which allows the calculation of a multiplicative factor 
connecting the observed luminosity to the SFR (for a dis- 
cussion on the conversion factors used in the literature at 
various wavelengths, see Kennicutt 1998a). The UV lumi- 
nosity density is dominated by the light emitted by young 
and massive stars, thus it provides an estimate of the cos- 
mic SFR density. Relevant observations of galaxies in the 
UV are typically performed at 1500 A ( Madau et al. 1996, 
Madau et al. 1998, Steidel et al. 1999, Massarotti et al. 
2001), 2000 A (Treyer et al. 1998) and 2800 A (Lilly et 
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al. 1996, Sawicki ct al. 1997, Comiolly ct al. 1997, Cowie 
et al. 1999). Other ways to assess star formation include 
observations of nebular emission lines such as Ha (Gallego 
et al. 1995, Gronwah 1998, Tresse & Maddox 1998, Glaze- 
brok ct al. 1999) and Oil (Hammer ct al. 1997). In all 
of these cases the main source of uncertainty in such de- 
terminations is still represented by interstellar extinction, 
caused either by the interstellar medium (ISM) of the ob- 
served galaxy or by our Galaxy. Star formation activity is 
likely to take place in highly obscured regions, so reliable 
extinction corrections are required to have accurate esti- 
mates of the SFR density based on observations in the UV 
band. 

However, without applying any extinction correction, the 
possibility that very intense star forming regions arc heav- 
ily obscured by dust implies that observations at short 
wavelengths can only provide lower limits to the actual 
SFR density. For this reason, the rc-processcd mid- and 
FIR- emission from the dust grains, which are heated by 
the UV and optical light emitted by the young stars, is a 
more reliable indicator of star formation activity, and it 
can be used to estimate the dust amount and SFR with 
fewer selection effects (Blain et al. 1999, Chary & Elbaz 
2001). The FIR studies can however be complicated by 
the multi-component effects of the dust (warm dust, cir- 
rus, Schaerer 1999) and by the old stars and AGN con- 
tributions to the heating of the dust. Observations in 
the FIR and submm bands include the results by Rowan- 
Robinson et al. (1997), Hughes et al. (1998) and Flores et 
al. (1999). More recently, observations in the radio con- 
tinuum have provided other estimates of the SFR at high 
redshift (Mobashcr ct al. 1999, Haarsma et al. 2000). Ob- 
viously, the discrepancies seen in the observed evolution of 
the LD are reflected in the measures of the SFR density, 
with a globally accepted increasing trend between z = 
and z 1, and a flat behaviour at higher redshifts, though 
the existence of a peak somewhere at 2 > 1 still cannot be 
completely excluded. 

3. CHEMICAL AND PHOTOMETRIC EVOLUTION OF 
GALAXIES 

By means of chemo-spectrophotometric models of 
galaxy evolution we aim at reconstructing the history of 
the luminous matter in the universe. These models al- 
low us to follow in detail the evolution of the abundances 
of several chemical species, starting from the matter re- 
processed by the stars and restored into the ISM through 
stellar winds and supernova explosions. We differentiate 
the galaxy types into ellipticals, spiral disks and irregulars. 
Wc assume that the category of galactic bulges is naturally 
included in the one of elliptical galaxies. Our assumption 
is motivated by the fact that they have very similar fea- 
tures: for instance, both are dominated by old stellar pop- 
ulations and respect the same fundamental plane (Binney 
& Merrifield 1998). This certainly indicates that they are 
likely to have a common origin, i.e. both are likely to have 
formed on very short timescales and a long time ago. 
Detailed descriptions of the chemical evolution models can 
be found in Mattcucci & Tornambe (1987) and Matteucci 
(1994) for elliptical galaxies, Chiappini et al. (1997) and 
Chiappini et al. (2001) for spirals and Bradamante et al. 
(1998) for irregular galaxies. It is worth noting that we 



choose these models with their specific assumptions be- 
cause they reproduce at best the observed local properties 
of the various galaxy types. 

According to our scheme, elliptical galaxies form as the 
result of the rapid collapse of a homogeneous sphere of 
primordial gas where star formation is taking place at the 
same time as the collapse proceeds, and evolve as "closed- 
boxes" , i.e. without any interaction with the external envi- 
ronment. Star formation is assumed to halt as the energy 
of the ISM, heated by stellar winds and SN explosions, 
balances the binding energy of the gas. At this time a 
galactic wind occurs, sweeping away almost all the resid- 
ual gas. 

Spiral galaxies are assumed to form as a result of two main 
infall episodes. During the first episode the halo forms and 
the gas shed by the halo rapidly gathers in the center yield- 
ing the formation of the bulge. During the second episode, 
a slower infall of external gas forms the disk with the gas 
accumulating faster in the inner than in the outer region 
("inside-out" scenario, Mattcucci & Francois 1989). The 
process of disk formation is much longer than the halo and 
bulge formation, with time scales varying from ~ 2Gyr in 
the inner disk to ^ 8Gyr in the solar region and up to 
15 — 20Gyr in the outer disk. 

Finally, irregular dwarf galaxies are assumed to assemble 
from merging of protogalactic small clouds of primordial 
composition, until a mass of ~ 6 x IO^Mq is accumu- 
lated, and to produce stars at a lower rate than spirals. 
In the next section, we will present a schematic outline of 
the equations and physical hypotheses at the basis of the 
models. 

3.1. The chemical evolution models 

Let Gi be the fractional mass of the element i in the 
gas within a galaxy, its temporal evolution is described by 
the basic equation: 

Gi = -4>{t)X,{t) + Ri{t) + {G,)raf - {G^)out (3) 

where Gi{t) = Mg{t)Xi{t) / Mtot is the gas mass in the 
form of an element i normalized to a total initial mass 
Mtot- The quantity Xi{t) = Gi{t)/G{t) represents the 
abundance in mass of an element z, with the summation 
over all elements in the gas mixture being equal to unity. 
The quantity G{t) = Mg{t)/Mtot is the total fractional 
mass of gas present in the galaxy at time t. ij){t) is the 
instantaneous SFR, namely the fractional amount of gas 
turning into stars per unit time; Ri{t) represents the re- 
turned fraction of matter in the form of an element i that 
the stars eject into the ISM through stellar winds and SN 
explosions; this term contains all the prescriptions regard- 
ing the stellar yields and the SN progenitor models. The 
two terms {Gi)inf and {Gi)out account for the inf ailing 
external gas from the intergalactic medium and for the 
outflow, occurring by means of SN driven galactic winds, 
respectively. The main feature characterizing a particular 
morphological galactic type is represented by the prescrip- 
tion adopted for the star formation history, summarized in 
the SFR expression. 

In the case of elliptical galaxies the SFR tp{t) (in Gyr~^) 
has a simple form and is given by: 

m = uG{t) (4) 
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The quantity v is the efficiency of star formation, namely 
the inverse of the typical time scale for star formation. In 
the case of ellipticals, v is assumed to drop to zero at the 
onset of a galactic wind, which develops as the thermal 
energy of the gas heated by supernova explosions exceeds 
the binding energy of the gas (Arimoto and Yoshii 1987, 
Matteucci and Tornambe 1987). This quantity is strongly 
influenced by assumptions concerning the presence and 
distribution of dark matter (Matteucci 1992); for the 
model adopted here a diffuse (i?e/-Rd=0.1, where iig is 
the effective radius of the galaxy and is the radius of 
the dark matter core) but massive {M^arklMLum = 10) 
dark halo has been assumed. 

In the case of irregular galaxies we have assumed a contin- 
uous star formation rate always expressed as in (4), but 
characterized by an efficiency lower than the one adopted 
for ellipticals. 

In the case of spiral galaxies, the SFR expression (Chiap- 
pini et al. 1997) is: 

Hr,t) = -[^^.]^('=-^^[^^^]'=-V,„.(r,t) (5) 
cr(r0,i) a[r,t) 

where i' is the SF efficiency, (j(r, t) is the total (gas + stars) 
surface mass density at a radius r and time t, cr(r0 , t) is the 
total surface mass density in the solar region and a gas {r, t) 
is the surface gas density. For the gas density exponent k a 
value of 1.5 has been assumed by Chiappini et al. (1997) in 
order to ensure a good fit to the observational constraints 
at the solar vicinity and in agreement with the estimates 
by Kennicutt (1998b). The three different star formation 
rates for ellipticals, spirals and irregulars as functions of 
time and for a Salpeter IMF are shown in figure 1. 

3.1.1. The Supernova Rates 

The true nature of Type la SNc and their progenitors 
is currently a matter of debate. Several models have been 
proposed so far, but divided between two main competing 
scenarios. According to the single degenerate (SD, Whe- 
lan & Iben 1973) scenario, a C-0 white dwarf accretes 
mass from a non-dcgencratc companion until it reaches 
the Chandrasekhar mass (~ 1.4M©) and explodes via C- 
deflagration leaving no remnant. The alternative double- 
degenerate (DD, Iben & Tutukov 1984) model sees the 
merging of two C-0 white dwarfs which, owing to the loss 
of angular momentum via gravitational wave radiation, co- 
alesce and explode by C-deflagration. For our purposes we 
adopt the SD model, which generally best represents the 
characteristics of the majority of type la SNe, as discussed 
in Matteucci & Recchi (2001). Type II SNe explosions are 
supposed to originate from core collapse of single mas- 
sive (M > 8M0) stars, with maximum masses allowed of 
IOOM0. 

The type la SN rate is expressed as: 

Rla{t) = A <^(Mb)[ / /(M)V'(i - TMMdMs 

where A is the fraction of binary systems which can end as 
SNIa within the IMF mass range, and which have a total 
mass AlBm < Mb < Mbm- A* = M2/MB is the ratio be- 
tween the secondary component of the binary system (i.e. 
the less massive one) and the total mass of the system and 



/(/x) is the distribution function of this ratio. Statistical 
studies indicate that mass ratios close to 0.5 are preferred, 
so the formula: 

/(m) = 2^+^(1 + 7)M^ (7) 

is commonly adopted (Matteucci & Recchi 2001), with 
7 = 2 as a parameter, tm-2 is the lifetime of the less mas- 
sive star in the system, which lives longer and therefore 
determines the time-scale for the explosion. The assumed 
value of A is ~ 0.1. This value assures that the present 
time SN rates are reproduced in each galactic type. For 
the masses M^m and Mbm we chose the values SM© and 
I6M0, respectively (sec Matteucci & Greggio 1986). 
The type II SN rate is simply expressed as: 

Rnit) = {l-A) J^l;"" </.(M)V(t - TM)dM+ 

/•lOO 

/ cj>{M)ij{t - TM)dM (8) 

J Mbm 

The IMF has been assumed to be constant in space and 
time and among galaxies. We have compared results cal- 
culated with two different IMFs, i. e. the Salpeter (1955) 
and the Scalo (1986). 



3.2. The spectro-photometric model 

The spectro-photometric calculations have been per- 
formed by means of the model by Jimenez et al. (1998). 
This model is based on the stellar isochrones computed by 
Jimenez et al. (1998) and the stellar atmospheric models 
by Kurucz (1992). The main advantage of this photo- 
metric code is that it allows one to follow in detail at 
every single timestep the mctallicity evolution of the gas 
out of which the stars form, at variance with other pop- 
ular photometric models which require the assumption 
of a constant metallicity. Using the stellar inputs, first 
we build simple stellar population (SSP) models consis- 
tent with the chemical evolution at any given time and 
weighted according to the assumed IMF. Then, a compos- 
ite stellar population (CSP) consists of the sum of different 
SSP formed at different times, with a luminosity at an age 
to and at a particular wavelength A given by: 



Lx{to) = r ( ' ^pit- to)Lssp,\{Z, t - to)dZdt (9) 

JO JZi 

where the luminosity of the SSP can be written as: 



LsspAZ,to-t) = / <t>{m)lx{Z,MM-t)dM (10) 

JMmin 

and 1\{Z, M, to — t) is the luminosity of a star of mass M, 
metallicity Z and age to — t; Zi and Zf are the initial and 
final metallicities, Mmin and M^ax are the smallest and 
largest stellar mass in the population, (j){m) is the IMF 
and ijj{t) is the SFR at the time t. 
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Fig. 1. — Star formation rates as functions of time for a spiral, elliptical and irregular galaxy in the case of a Salpeter IMF calculated 
according to equations (4) and (5). 



4. RESULTS 

4.1. The evolution of the galaxy luminosity density 

In the B band, at ^; = the LDs for the single galaxy 

types are simply given by the integral of the LFs observed 
by Marzke et al. (1998, see their equation 2). At redshift 
other than zero we tranform the absolute magnitudes ap)- 
plying the evolutionary corrections, calculated by means 
of the spectro-photometric code for every galaxy type: 



J Ex/i+Az)RB(\)d\ 

'''"'"''-°'+^-^"'»' }^v:::(oU(aU ' 



(11) 

where Mb{z = 0) and Mb{z) are the absolute blue mag- 
nitudes at redshift and z, respectively, Ex{z)dX is the 
energy per unit time radiated at the rest-frame wavelength 
A by the galaxy at redshift z, and Rb{X) is the response 
function of the rest-frame B band. The second term on 
the right side of equation (11) represents the evolutionary 
correction (EC), i.e. the difference in absolute magnitude 
measured in the rest frame of the galaxy at the wavelength 
of emission (Poggianti 1997). 

We then calculate the B band LF at redshift z according to: 



^b{Mb,z) = ^b{Mb{z)) 



(12) 



which is equivalent to assuming evolution in luminosity 
and not in number. This clearly coincides to assume that 
the effects of mergers are small at all redshifts. In bands 
other than B we assume that the LF shape is the same 
as in the B band and we calculate the LF in the given 
band (X) transforming the absolute magnitudes according 
to the rest-frame galaxy colors as predicted by the spec- 
trophotometric model: 



Finally, the total luminosity density per unit frequency in 
a given band px is given by the sum of the single contri- 
butions of the three different galactic morphological types 
(ellipticals, spirals and irregulars): 



(14) 



Mx=Mb + {X- B)rf 



(13) 



Figures 2 and 3 show the predicted evolution of the to- 
tal LD in the U (centered at 3650A), B (centered at 4450 
A), I (centered at 8060 A) and K (centered at 21900 A) 
bands and the contributions of each morphological type for 
a Salpeter (1955) and a Scalo (1986) IMF, respectively. In 
all bands, at early times the total LD is dominated by 
the light produced by elliptical galaxies, which experience 
their strong star-burst during which they are the most in- 
tense sources in each photometric band. After this phase, 
which lasts ~ 0.3 Gyr, they will evolve passively for the 
rest of their life. On the other hand, the disks of spiral 
galaxies form stars continuously at all epochs: their lumi- 
nosity increases slowly until 3Gyr from the beginning of 
the star formation, where it reaches a broad peak, then 
it starts to decrease significantly, owing to the progressive 
consumption of their gas reservoirs. 

Irregular galaxies are the most slowly evolving systems and 
their star formation rates never reach the values recorded 
in spiral disks and ellipticals. As a consequence, their LD 
values are the smallest at any time and in any band. 

In the I and K bands the difference between the spiral 
and elliptical luminosity densities are small after 3 Gyr. 
At this epoch the total amount of stars is nearly in place 
in both types. At the present time (13 Gyr), the con- 
tributions of ellipticals and spirals to the total luminos- 
ity density in the I and K bands is practically the same, 
at variance with what happens in the B and U bands. 
The difference between the Salpeter and the Scalo IMF 
concerns the stars in the high mass regime (M > 8Mq), 
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Fig. 2. — Predicted evolution of the luminosity density for galaxies of different types in the U, B, I and K bands for a Salpeter IMF. Dotted 
line: elliptical galaxies; short-dashed line: spiral galaxies; long-dashed line: irregular galaxies. 
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Fig. 3. — Predicted evolution of the luminosity density for galaxies of different types in the U, B, I and K bands for a Scalo IMF. Dotted 
line: elliptical galaxies; short-dashed line: spiral galaxies; long-dashed line; irregular galaxies. 
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which are more numerous in the former case. This trans- 
lates into the amount of UV and B light emitted by all 
galaxies, which is sensitively larger with the Salpeter than 
with the Scalo IMF. At longer wavelengths the difference 
diminishes: in the I band, where there is still a significant 
contribution brought by high intermediate mass stars ex- 
periencing the asymptotic giant branch (AGB) phase, the 
difference between the elliptical and spiral luminosity is 
stronger with the Scalo IMF, whereas the K band lumi- 
nosities have practically the same behaviour in both c;ases. 

Figure 4 shows the predicted and observed evolution of 
the LD in various bands for a Salpeter (left panels) and 
a Scalo (right panels) IMF and for a galaxy formation 
redshift of Zf = 5, having assumed an Einstein-De Sitter 
(EdS) cosmological model (Oq = 1, Oa = 0) and h = 0.5. 
Figure 5 shows the same as in figure 4 but in the case of 
a Lambda-cold dark matter cosmology (ACDM, Qq = 0.3, 
fl\ = 0.7) and h = 0.65. Figures 6 and 7 show the same 
as in figures 4 and 5, respectively, but with a galaxy for- 
mation redshift of z/ = 10. The data have been converted 
from the EdS to the ACDM cosmology according to the 
prescriptions given by Somerville et al. (2001). Given the 
slight discrepancy between the B-band LF normalization 
by Marzke et al. (1998), adopted for our estimate in the 
local universe, and the one by Ellis et al. (1996), which 
is shown in the plots, we have re-normalized our value at 
z = to the one by Ellis et al. (1996) in order to better 
stress the goodness of our fit. 

In the U band (upper panels) the data show a consid- 
erable spread, with large error bars at any redshift. All 
the UV data reported in the figure are not extinction- 
corrected, with the only exception of the three points mea- 
sured by Massarotti et al. (2001). Our predictions indicate 
a peak in the total UV LD corresponding to the redshift 
of galaxy formation and a monotonically decreasing trend 
down to z = 0, with a particularly strong evolution be- 
tween z — 1 and z = 0. This last feature is in agreement 
with the data by Pascarelle et al. (1998) and with the 
low redshift value by Treyer et al. (1998). We do not see 
evidence for the very steep drop observed by Lilly et al. 
(1996). From figures 4 and 5 we note that, if the peak 
generated by elliptical formation were shifted at z ~ 3, 
it would be fairly consistent with the extinction-corrected 
measures by Massarotti et al. (2001), in particular in the 
case of the Scalo IMF. To our knowledge, the extinction- 
corrected data by Massarotti et al. (2001) represent the 
first observational evidence for a very high peak in the 
observed UV LD which could be clearly ascribed to the 
massive star formation occurring during the formation of 
spheroids. 

In the B-band (lower panels) there is a smaller number of 
detections than in the UV but a better agreement between 
the observations and the predictions in the < z < 1 red- 
shift range. The smallest number of observations is found 
in the I and K bands. At z = the K-band LD is rather 
well constrained by several measures (Gardner et al. 1997, 
Loveday 2000, Cole et al. 2001, Kochanek et al. 2001), 
which show a substantial agreement with one another. For 
clarity, in our plots we have chosen one single value, the 
one by Gardner et al. (1997). At redshift larger than 
zero, in the I band the only observations at our disposal 
are the ones by Lilly et al. (1996) at 1 /x. We note that 



our curves reproduce very well the oserved evolution of the 
LD in the B, I and K bands. The best agreement between 
the data and the predictions is achieved in the case of the 
ACDM cosmology, especially in the I and K bands. This 
occurs mainly because, having assumed galaxy formation 
redshifts of z/ = 5 and z/ = 10, in the ACDM cosmology 
a formation redshift of z/ =5 (z/ = 10) corresponds to a 
lookback time of 13.26 Gyr (14 Gyr), whereas in the EDS 
Zf = 5 {zf = 10) corresponds to 12.15 Gyr (12.68 Gyr). 
This means that in the ACDM galaxies are older than in 
the EDS, consequently they contain redder stellar popu- 
lations. Both the Salpeter and the Scalo IMFs provide 
satisfactory results. As already mentioned, the Salpeter 
IMF is richer in high mass stars than the Scalo IMF. The 
consequence of this is that with the Salpeter IMF all galax- 
ies are more luminous at the shortest wavelengths, where 
these stars emit the bulk of their light. Moreover, with the 
Salpeter IMF galaxies are on average redder than with the 
Scalo one. This translates into values of the I and K LDs 
slightly higher. However, given the uncertainties and the 
small data sample in the I and K bands, it is difficult to 
conclude which choice for the IMF should be preferred. 
Our main result is that the disks of spiral galaxies are 
mainly responsible for the observed decrease of the total LD 
in the U and B bands between z ^ 1 and z ~ 0, whereas 
in the I band the decline is due to contributions from both 
spirals and ellipticals, as can be seen also in figures 2 and 
3. 

Figures 6 and 7 show the evolution of the comoving LD 
in various bands assuming that galaxies formed at z/ = 10, 
having assumed an EdS and ACDM cosmology, respec- 
tively, and in the cases of a Salpeter (left panels) and 
a Scalo (right panels) IMF. The predictions are still in 
fair agreement with the observations, in particular in the 
< z < 1 redshift range. The main differences concern 
the predictions at high redshift: the Zf = 10 scenarios 
allow for a nearly constant LD throughout a wider red- 
shift range. As far as the z < 1 behaviour of the LD is 
concerned, little difference can be seen when we assume 
a global galaxy formation atz/ = 5orz/ = 10. This is 
due to the short cosmic time elapsed between z = 5 and 
z = 10 both in the EdS and ACDM cosmologies, namely 
0.53 and 0.74 Gyr, respectively. This time scale is much 
shorter than the one corresponding to the < z < 1 red- 
shift range, i. e. 8.43 Gyr in the EdS and 8.31 Gyr in 
the ACDM, having thus little influence on what happens 
during the latter period. 

The study of the galaxy LD at high redshift is crucial for 

understanding the epoch of spheroid formation. Unfortu- 
nately, the high redshift measures concern only the UV 
band, which is the most seriously affected by dust extinc- 
tion effects (Adelberger & Steidel 2000). If the attenua- 
tion by dust is as serious as advocated by Massarotti et 
al. (2001), the peak generated by the formation of mas- 
sive spheroids could lie anywhere in the z > 1 region, but 
it would be hidden if the associated strong starbursts oc- 
curred in sites heavily obscured by dust. In fact, spheroids 
are the systems which most rapidly reach oversolar metal- 
hcities (Pettini et al. 2002, Matteucci & Pipino 2002), 
therefore the formation of dust grains could be particu- 
larly favoured if, as is reasonable to assume, the proba- 
bility of dust formation is proportional to the metal con- 
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Fig. 4. — Upper Panels: predicted (solid line) and observed luminosity densities in the U band for a Salpeter (left panel) and Scalo (right 
panel) IMF in the case of galaxy formation at Zf = 5. Lower panels: predicted and observed luminosity densities in the B, I, K bands. The 
adopted cosmology is the one of Einstein-De Sitter (SIq = 1, Oj\ = 0) with h = 0.5. Data in the upper panels: Solid hexagons from Pascarelle 
et al. (1999, 1500 A), open hexagons from Steidel et al. 1999, cross from Treyer et al. (1998, 2000 A), stars from Massarotti et al. (2001, 
1500 A), solid pentagons from Madau et al. (1998, 1500 A), solid circles from Lilly et al. (1996, 2800 A), open squares from Connolly ct al. 
(1997, 2800 A, four-points stars from Cowie et al. 1999. The three-points stars are the extinction-corrected data by Massarotti et al. (2001). 
Data in the lower panels: Open triangles from Ellis et al. (1996, 4400 A, B), open pentagons from Lilly et al. (1996, 4400 A, B), solid squares 
from Connolly et al. (1997, 4400 A, B), open circles from Lilly et al. (1996, 1 ^j, I), solid triangle from Gardner et al. (1997, 2.2 ^, K). 




Fig. 5. — Upper Panels: predicted (solid line) and observed luminosity densities in the U band for a Salpeter (left panel) and Scalo (right 
panel) IMF in the case of galaxy formation at = 5. Lower panels: predicted and observed luminosity densities in the B, I, K bands. The 
adopted cosmology is the ACDM (Cq = 0.3, Ha = 0.7) with h = 0.65. Data are as described in figure 4. 
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Fig. 6. — Upper Panels: predicted (solid line) and observed luminosity densities in the U band for a Salpeter (left panel) and Scalo (right 
panel) IMF in the case of galaxy formation at Zf = 10. Lower panels: predicted and observed luminosity densities in the B, I, K bands. The 
adopted cosmology is the one of Einstein-De Sitter {Uq = 1, Ha = 0) with h = 0.5. Data are as described in figure 4. 




Fig. 7. — Upper Panels: predicted (solid line) and observed luminosity densities in the U band for a Salpeter (left panel) and Scalo (right 
panel) IMF in the case of galaxy formation at Zf = 10. Lower panels: predicted and observed luminosity densities in the B, I, K bands. The 
adopted cosmology is the ACDM (Co = 0.3, Ha = 0.7) with h = 0.65. Data are as described in figure 4. 
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tent. High redshift observations at higher wavelengths, 
i. e. in bands not affected by biases caused by dust ex- 
tinction, could provide fundamental hints on the epoch of 
major spheroid formation, as well as high redshift obser- 
vations in the FIR/sub-mm bands, where dust rc-cmits all 
the starlight absorbed in the rest-frame B and U bands 
(Hughes et al. 1998, Blain et al. 1999). 

4.2. The cosmic star formation rate density 

Figure 8 shows the evolution of the cosmic SFR density 
as a function of redshift (Madau's plot) as predicted by 
our models and as observed by several authors in the case 
of a ACDM cosmology and galaxy formation at Zf = 5. 
The cosmic SFR density is not a directly observed quan- 
tity: to be evaluated it needs, besides the measure of the 
LD at certain wavelengths (mostly 1500 and 2000A), the 
adoption of a universal IMF to calculate the conversion 
factor between the observed p\ and p*. For this reason 
the determination of this quantity is affected by two main 
sources of uncertainty. One is related to the determination 
of the calibration constant, uncertain if the universal IMF 
is varying throughout cosmic time. The second is related 
to the dust extinction affecting the UV observations. We 
have computed the cosmic SFR density according to: 

Mz) =^PBiiz){^)B^iz)^p^iz) (15) 
i 

where pBi, (x)si ^^'^ V'i ^-re the B LD, the B mass-to- light 
ratio and star formation rate for the galaxies of the i-th 
morphological type, respectively. The points in figure 8 
represent estimates based on observations at short wave- 
lengths, except the solid hexagon obtained by means of 
sub-mm measures by Hughes et al. (1998). However, since 
the extent of the attenuation by dust is highly uncertain 
(see Steidel et al. 1999, Hopkins et al. 2001, Thompson 
et al. 2001), we have chosen to use the data uncorrected 
for dust extinction. For this reason, all the data based 
on optical observations should be regarded as lower limits 
to the true values. The thick lines in the curves indicate 
lower and upper limits on the SFR density calculated by 
Haarsma et al. (2000) from radio observations. As well 
as the point from sub-mm data, these estimates need no 
corrections for dust obscuration since the radio emission 
at V > \GHz passes freely through dust. 
At 2 ~ we overpredict the value observed by Gallego 
et al. (1995) by a factor of 4. Interestingly, our estimate 
is instead in very good agreement with a recent determi- 
nation based on Local Group observations by Hopkins, 
Irwin and Connolly (2001), which calculate a local SFR 
density higher than the one by Gallego et al. (1995) by 
a factor of ~ 5. This fact implies that, even in the local 
universe, dust obscuration could introduce serious biases 
in the determinations of the global SFR density. Since the 
probability to intercept dusty objects increases in lockstep 
with the line of sight, it is conceivable that this effect could 
become more and more pronounced at increasing redshift. 
Another source of discrepancy between our local value and 
the one by Gallego et al. (1995) is the normalization of the 
luminosity density that we have chosen, which in this case 
seems too high. Since the discrepancy between our local 
value and the observed one arises from the combination 



of both effects, it seems very difficult to quantify the un- 
certainty associated to the choice of the normalization. If 
we were to neglect all the effects due to dust obscuration, 
and if we took the value by Gallego et al. (1995) at face 
value, the uncertainty due to the normalization would be 
of a factor of 4. In figure 8, in order to better stress the 
comparison between our prediction and the observations, 
our curve has been re-normalized to the value by Gallego 
et al. (1995). Our predictions indicate a peak in the SFR 
density at high redshift due to massive star formation in 
spheroids. In our models, the SFR in ellipticals can reach 
very high values (up to lOOOMo/yr), very similar to the 
ones observed in some SCUBA (Ivison et al. 2000) and 
Hyperlouminous infrared galaxies (Rowan-Robinson et al. 
1997, Rowan-Robinson 2000). Such values are motivated 
by the fact that, in order to reproduce several observa- 
tional features such as the observed increase of the [Mg/Fe] 
ratio and the velocity dispersion cr, giant elliptical galax- 
ies, i.e. with masses of the order of ^ IQ^^Mq, have to 
form all their stars on very short timescales, i.e. r < 1 
Gyr (Matteucci 1994). Hence, the possible SFR rates for 
high mass ellipticals span the range ~ 100 — lOOOMo/t/r. 
If we were to assume a galaxy of M ^ 10^°Mq as a typical 
elliptical, with a star formation timescale of ~ 1 Gyr, this 
would lower the peak in the SFR density by a factor of 
10. As already mentioned in section 4.1, if most of the 
star formation in the universe occurred at very high red- 
shift and in sites highly obscured by dust, there could be a 
peak somewhere &t z > 1 which would remain completely 
unseen in the observed />*-z plot. 

4.3. The evolution of the stellar mass density 

Figure 9 shows the predicted and observed evolution of 
f2* , namely the stellar mass density (living stars plus rem- 
nants) divided by the critical density of the universe. The 
stellar mass density is given by: 

P*{z) = ^pBi{z){—)Bi{z) (16) 

i 

where PBi{z) is the B luminosity density and (^)Bi is the 
B mass-to-light ratio for the galaxy of the i—th morpho- 
logical type. In figure 9 we show two values of the stellar 
mass density measured in the local universe: the one by 
Cole et al. (2001), valid for the whole galaxy population, 
and the one by Perez-Gonzalez et al. (2003), repesenting 
the contribution brought by local star forming galaxies. 
The discrepancy between our results and the observations 
(and among the observations) at 2 = is likely to be as- 
cribed to uncertainties in the normalization chosen for the 
local LF. We notice that our prediction for the total stellar 
mass density is consistent with the values by Brinchmann 
& Ellis (2000) and with the low redshift value by Dickin- 
son et al. (2003). In our scenario, ellipticals build up the 
totality of their stellar mass at very early times, whereas 
disks of spirals and irregulars build up their stars progres- 
sively. Thus, the rise in the 3 > z > 0.5 redshift range, 
observed by Dickinson et al. (2003), is mainly determined 
by spiral galaxies, namely by the same galaxies contribut- 
ing to the drop in the cosmic SFR density between z ^ 1 
and z = 0. 

Finally, it is worth noting that our z = mass density val- 
ues are in good agreement with the estimates by Fukugita 
et al. (1998) for spheroids, disks and irregulars. 
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Fig. 8. — Predicted and observed star formation rate density (left scale) and metal-ejection rate density (right scale) for a Salpeter IMF 
in the case of a ACDM cosmology and galaxy formation at Zf = 5. All the points are taken from Somerville et al. 2001. The thick lines 
represent lower and upper limits based on radio observations (Haarsma ct al. 2000). 




Fig. 9. — Predicted and observed comoving stellar mass density f2* for a ACDM cosmology, a Salpeter IMF and galaxy formation occurred 
at Zf = 5. Dotted line: predicted evolution of the stellar mass density (living stars plus remnants) in elliptical galaxies; short-dashed line: 
predicted contribution from spiral galaxies; long-dashed line: predicted contribution from irregular galaxies; solid line: predicted total evolu- 
tion of the stellar mass density. The boxes represent the values observed by Dickinson et al. (2003), the open circle is the 2 = value observed 
by Cole et al. (2001), the open triangles are by Brinchmann & Ellis (2000) and the open square is the value observed by Perez-Gonzalez et 
al. (2003) in local star-forming galaxies. The solid triangle, solid square and solid pentagon are the local estimates by Fukugita et al. (1998) 
of the stellar mass density in spheroids, disks and irregulars, respectively. 
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4.4. The missing metals crisis 

The metals in the high redshift universe are observed 
in various amounts in different sites. One is represented 
by the DLA systems, which are considered to be the high- 
redshift counterparts of the gas-rich galaxies in the local 
universe and are usually associated with slowly-evolving 
galactic or sub-galactic systems (Pettini et al. 1999, 
Prochaska et al. 1999, Calura et al. 2002). Other sites 
of significant metal production at high redshift are the 
Lyman-break galaxies (LBG, Steidel et al. 1996, Steidel 
et al. 1999) whose high metal abundances, strong lumi- 
nosities and kinematical features indicate their possible as- 
sociation with star-forming spheroids (Matteucci & Pipino 
2002, Pettini et al. 2002). A certain amount of the metals 
produced in proto-galaxies is ejected through SN-driven 
winds into the inter-galactic medium, whose presence is 
detected as a "forest" of absorption lines bluewards of the 
Lyman-a emission lines of high redshift QSOs (for a re- 
view of the argument, see Pettini 2000, Bechtold 2001). 
The inter-galactic gas is highly ionised and is likely to ac- 
count for most of the baryons at low and high redshift. In 
this medium the lowest metallicity in the universe can be 
observed (Songaila 2001), possibly due to enrichment by 
pop III stars (Ostriker & Gnedin 1996). 
Finally, an uncertain amount of metals is locked up in 
QSO and active galactic nuclei (Hamann 1997), objects 
with metallicities even higher than the ones observed in 
LBGs (Matteucci & Padovani 1993). 
We evaluate the metal ejection rate (MER) density ac- 
cording to Pagel (2001) by assuming the proportionality 
between the comoving density of metals produced by stars 
in the universe and the SFR density to be a factor 1/42 
(Madau et al. 1996, Pettini 1999, Pagel 2001). This is 
in principle a rough approximation, since it is equivalent 
to ignore the contribution of type la SNe and consider 
only the elements produced explosively by massive stars. 
However, for our purposes this assumption is reasonable 
since the main contributor to the total metallicity Z is the 
oxygen, which is produced on very short timeseales exclu- 
sively by massive stars. A more refined calculation of the 
production rate in the universe of various elements will be 
the object of a forthcoming paper (Calura & Matteucci, 
in preparation). 

If we integrate the predicted MER dc;nsity across the 
epoch ranging from 13.3 and 11 Gyr, which corresponds, 
for the ACDM cosmology, to the redshift range where the 
bulk of DLAs and LBG are found, we obtain: 



13.3Gyr 



llGyr 



pz{t)dt = 5.2 X IO^MqMpc- 



(17) 



An approximative estimate based on the observed values of 
p* yields 5 x 10^ MqMpc~^ (Pagel 2001), a value in excel- 
lent agreement with our appraisal, though slightly lower. 
The observed contribution from DLAs and LBGs repre- 
sents ~ 10% of the estimate by Pagel and us. We find 
such a high value for the total metal density at high red- 
shift since we assume that all ellipticals form at the same 
time. If we spread the formation of ellipticals over a wide 
redshift range, we would probably find a lower value both 
for the SFR density peak and for the metal density. For 
this reason, the values predicted for such peak and for the 



total amount of missing metals at high redshift should be 
regarded as upper limits. However, this does not mean 
that the missing metal crisis could not be even more se- 
rious than what the observations indicate if most of the 
metals were produced in dust-obscured starbursts asso- 
ciated with early spheroid formation. We suggest that, 
irrespective of the amount of metals, the most plausible 
site of the missing metals could be the warm gas in galaxy 
groups and proto-clusters, in which the metals could have 
been ejected through strong winds following the intense 
starbursts in ellipticals. In such an environment, the pres- 
ence of diffuse metals could be difficult to detect, likely 
owing to virial temperatures lower than the ones observed 
in clusters (Renzini 1997, Fukugita et al. 1998). The main 
sources of strong winds could be the LBGs, where recent 
kinematic studies have indicated the frequent presence of 
strong large-scale outfiows (Pettini et al. 2001, 2002). An- 
other possible source is represented by the SCUBA (Blain 
et al. 1999, Trentham, Blain & Goldader 1999) and hyper- 
luminous infrared galaxies (Rowan-Robinson 2000), very 
luminous in the infrared band but faint in the optical, 
where the most intense observed starbursts are likely to 
occur and which are natural candidates for massive proto- 
ellipticals (Ivison et al. 2000). 

4.5. The cosmic SN rate 

Another key-issue in the study of galaxy evolution is the 
determination of the cosmic SN rate, which can provide 
useful constraints on the cosmological parameters (Marri, 
Ferrara fc Pozzetti 2000), on the possible SN progenitor 
models and on the evolution of the cosmic star formation 
history. The observed cosmic SN rate is expressed in SNu, 
i. e. in SN per unit time per unit blue luminosity: 



ISNu = ISN/WLqb/ century 



(18) 



The observed SN rate R is expressed as (Hardin et al. 
2000): 

R = N/S (19) 

where N is the number of SN detected and S is a sum over 
the observed galaxies weighted by their blue luminosities 
Li-. 



ei{t,Zi)dt 



(20) 



where e is the efficiency to detect in the i—th galaxy a SN 
whose maximum occurs at time t in the SN rest frame. 
For similarity, we calculate the cosmic SN rate at redshift 
z as: 

J2^ ri{z) 



Rc{z) 



(21) 



where ri{z) represents the number of SNe per 100 yr ex- 
ploding in the i—th galaxy type at redshift z, whereas Lbi 
is the predicted blue luminosity of the i—th galaxy type 
at redshift z. Figures 10 and 11 show the results of our 
computation of the cosmic type la and II SN rates, cal- 
culated in SNu (1 SNu= 1 SN per century per IO^^Lqb), 
having assumed a ACDM world model and galaxy forma- 
tion occurring at = 5. In figure 10 we adopt a universal 
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Fig. 10. — Predicted and observed type la (below) and II (above) SN rate expressed in SNu (SN ccntury^^ 10^^" Lqb) for a Salpeter 
IMF, a ACDM cosmology and galaxy formation at Zf = 5. Dotted line: elliptical galaxies; short-dashed line: spiral galaxies; long-dashed 
line: irregular galaxies; solid line; total SN rate. Open squares: Cappellaro et al. 1999; open pentagon: Hardin et al. 2000; solid pentagon: 
Pain et al. 1996; solid square: Pain et al. 2002. 




redshift 



Fig. 11. — Predicted and observed type la (below) and II (above) SN rate expressed in SNu (SN century ^ 10 ^'^Lqb) for a Scalo IMF, a 
ACDM cosmology and galaxy formation at = 5. Curves and symbols are as indicated in figure 10. 
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Salpctcr IMF, whereas in figure 11a Scalo IMF. We com- 
pare our predictions for galaxies of different morphologies 
with data observed by several authors up to 2; ~ 0.55 (see 
the captions to figures 10 and 11 for details). According 
to our calculations, in the local universe the bulk of SNc 
la explode in elliptical and spiral galaxies. At increasing 
redshift the contribution of ellipticals becomes more and 
more significant. At 2; ~ 0, practically all type II SNe ex- 
plode in spiral galaxies, with a small contribution brought 
by irregular galaxies. In ellipticals the type II SNe explode 
only at high redshift, since these galaxies are assumed to 
evolve passively since the interruption of the star forma- 
tion at the onset of the galactic wind, which occurs at 
early times. We note that our fit to the observed points 
is excellent in the case of the Salpeter IMF both for type 
la and II SNe. In the case of the Scalo IMF we under- 
estimate the number of SN II in the local universe and 
we slightly underpredict the type la SN rate observed at 
high redshift (z - 0.55, Pain et al. 2002). We stress that, 
unlike all our previous predictions, the calculation of the 
SN rate is independent on the luminosity function and is 
directly comparable to the observations without requiring 
a normalization at z = 0. For this reason, the calculation 
of the SN rate can be considered a very reliable test of 
consistency for our chemo-spectrophotometric model. 
If the cosmic SN rate can be considered as a reliable indica- 
tor of how star formation has evolved in the universe since 
the growth of galactic structures, as discussed in Madau, 
Delia Valle & Panagia (1998) and Sadat et al. (1998), our 
result is a further evidence that the bulk of the galaxy 
population has evolved in luminosity and not in number 
since early times, i.e. that the bulk of the galaxies as we 
see them today were already in place at early epochs. Ob- 
servations of the cosmic SN rate at very high redshifts, 
hopefully achievable after the launch of the Next Gener- 
ation Space Telescope (NGST), can provide fundamental 
clues on the galaxy formation epoch. 

5. CONCLUSIONS 

In this paper we have used detailed chemo- 
spectrophotometric models for galaxies of different mor- 
phological types, namely ellipticals, spirals and irregulars, 
in order to study the evolution of the galaxy luminosity 
density in various bands, the cosmic SFR density and 
the cosmic SN rate. We have studied the evolution of 
the baryonic matter in the universe in the form of stars 
and metals, starting from the B-band luminosity function 
observed in the local universe by Marzke et al. (1998), 
having assumed that the whole galaxy population started 
forming stars at high redshift. We have considered two 
different forms for the universal IMF, namely the Salpeter 
(1955) and the Scalo (1986), and we have computed our 
models for different cosmologies, i. e. the Einstein-De 
Sitter and the ACDM. We have compared our predictions 
with an updated large set of high-redshift observations 
and we found an overall good agreement between our pre- 
dictions and all the observational evidence considered in 
the present work. In particular, our main results are: 
1) We reproduce the evolution of the galaxy luminosity 
density in the B, I and K bands considering that all galaxy 
types evolve in luminosity and not in number since the 
epoch of their formation. Generally, we overestimate the 



luminosity density observed in the UV, in particular at 
high redshift. In our scheme, elliptical galaxies suffer 
high-redshift starbursts after which they evolve as passive 
systems, i. e. without star formation, whereas spiral disks 
and irregular galaxies form stars continuously down to the 
present epoch. At high (> 2) redshift, ellipticals dominate 
the total emissivity in any optical band. We predict a high 
redshift peak generated by the intense starburst in ellip- 
ticals which is not visible in the available observations at 
short wavelengths not corrected by extinction. The cause 
could be a heavy attenuation by interstellar dust which 
could hide the most intense star-forming sites, or the bulk 
of the rest-frame UV light emitted by single objects. On 
the other hand, extinction-corrected high redshift mea- 
sures of the UV flux density (Massarotti et al. 2001) 
are fully compatible with our predictions. High-redshift 
observations in bands less sensitive to dust attenuation, 
currently absent, would be fundamental in probing the 
existence of such a peak. At z < 2, in the U and B bands 
the disks of spiral galaxies emit the bulk of the light and 
are the main contributors to the decline in the galaxy 
luminosity density observed between 2; ~ 1 and z = 0, 
whereas in the the I and K bands the light emitted by 
stars in ellipticals is comparable to that emitted by the 
spirals. Irregular galaxies bring a negligible contribution 
to the total luminosity density at any epoch and in any 
band. Both the Salpeter and Scalo IMFs provide very 
satisfactory results, especially assuming a ACDM cosmo- 
logical model. 

2) We predict that spiral disks are the most active star 
forming sites in the universe throughout the redshift in- 
terval < z < 1 which, for a standard ACDM cosmology 
and h = 0.65, correspond to the 57% of the cosmic time. 
Furthermore, if the SFRs recorded in SCUBA and hyper- 
luminous infrared galaxies can be considered typical for 
starbursting protoellipticals, we predict a very high peak 
(up to ^ 1000 times the value observed in the local uni- 
verse) in the cosmic SFR density occurring at the epoch 
when the bulk of the spheroids formed. 

3) According to our calculations, the missing metal crisis 
coiild be even more serious than what has been proposed 
by previous authors (Pettini 1999, Pagel 2001) if the bulk 
of the metals were produced in dust-obscured starbursts 
associated with an early spheroid formation. We suggest 
that, regardless of the amount of metals, the most plau- 
sible site for the missing metals (;ould be the warm gas 
in galaxy groups and proto-clusters, in which the metals 
could have been ejected through strong winds following 
intense starbursts. In such an environment, the presence 
of the metals would be difficult to detect owing to the 
low virial temperatures. In principle, it could be possi- 
ble to detect these metals by identifing the most recently 
assembled clusters of galaxies, provided that they have 
formed out of groups where the missing metals lie. The 
main sources of such strong winds could be both Lyman- 
break galaxies (Pettini et al. 2002), SCUBA (Trentham 
et al. 1999) and hyperluminous infrared galaxies (Rowan- 
Robinson 2000), which are the best candidates for the 
high-redshift counterparts of the nearby massive elliptical 
galaxies. 

4) The study of the cosmic SN rate is a fundamental test 
of consistency for our models since it is completely inde- 
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pendent of the normalization of the galaxy population at 
z = 0, which was required by all our previous calculations. 
We reproduce very well the observed cosmic SN rate as- 
suming a universal Salpeter-like IMF. The assumption of 
a Scalo IMF causes an underestimate of the local type II 
SN rate and of the la SN rate observed at high redshift. 
We predict that most of the SN la explode in ellipticals 
at all ages, whereas all the SN II explode in spiral disks 
and irregulars after the halt of star formation in elliptical 
galaxies occurred at high redshift. This is in fair agree- 
ment with all our previous results, which indicate that 
spiral galaxies are the most active star forming objects 
since at least z ^ 1. Also the observed cosmic SN rate 
can be explained by assuming a scenario of coeval galax- 
ies whose number is conserved and which evolve purely 
in luminosity. The determination of the cosmic SN rate 
at redshifts > 0.55 (i.e. the highest explored by current 
SN surveys, Pain et al. 2002), hopefully achievable after 
the launch of the Next Generation Space Telescope, can 
provide further constraints on galaxy evolution models. 
All the results obtained in this paper clearly indicate that 
the bulk of galaxies was already in place at early epochs. 



According to several recent observations, the presence of 
high mass galaxies at high redshift is not rare (Saracco 
et al. 2002, Rusin et al. 2002, Pozzetti et al. 2003). 
Furthermore, in the last few years increasing evidence has 
been found that most (up to 70%) of the extremely red ob- 
jects discovered in deep NIR and optical surveys could be 
high-redshift ellipticals (Daddi et al. 2000 and references 
therein). Certainly, all these evidences seem difficult to 
reconcile with galaxy formation scenarios where the bulk 
of the stellar mass is formed at ^ < 1. 
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